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ff^ ' Abstract 

Supernovae searches have shown that a simple matter-dominated and decelerating universe should be ruled out. 
However a determination of the present deceleration parameter go through a simple kinematical description is not 
exempt of possible drawbacks. We show that, with a time dependent equation of state for the dark energy, a bias is 
present for go '■ models which are very far from the so-called Concordance Model can be accommodated by the data 
and a simple kinematical analysis can lead to wrong conclusions. We present a quantitative treatment of this bias 
and we present our conclusions when a possible dynamical dark energy is taken into account. 



PACS Numbers : 98.80.Es, 98.80. Cq 

Key- Words : cosmological parameters - supernovae 

Number of figures : 3 

February 2005 
CPT-2005/P.007 

anonymous ftp or gopher : cpt.univ-mrs.fr 
E-mail : virey@cpt.univ-mrs.fr 



* "Centre de Physique Theorique" is UMR 6207 - "Unite Mixte de Recherche" of CNRS and of the Universities "de 
Provence", "de la Mediterranee" and "du Sud Toulon- Var"- Laboratory affihated to FRUMAM (FR 2291). 

+ "Centre de Physique des Particules de Marseille" is UMR 6550 of CNRS/IN2P3 and of the University "de la 
Mediterranee" . 



1 



Observations 0, of Type la Supernovae (SNe) al- 
low to probe the expansion history of the universe. The 
measurements of the apparent magnitudes m{z) of SNe 
determine the luminosity distance versus redshift dL{z) 
which reads for a flat (SIt = 1) universe : 



dLiz) = c{l + z) 



du 



H{u) 

+ t\l.. A- S^[l+q{x)]d\n(l+x)\ (^^^ 



Ho 



du e I 



where dL{z) is written here in terms of the epoch- 
dependent deceleration parameter q(z) = {—a/ a) / H'^{z) 
with H{z) = a/ a. dL{z) is related to the measured mag- 
nitude by m{z) = Ms + 5logio{Ho/cdL{z)), Ms being 
a normalisation parameter which combines the absolute 
SNe magnitude and the Hubble constant. 

The present behavior of the expansion is attributed 
to a new "Dark Energy" (DE) component with negative 
pressure : px = w px , with a possibly time dependent 
equation of state (EoS) w{z) (for a recent review see Ref. 
13 )• fix (^m) will denote the ratio of the present DE 
(Matter) density to the critical density. In this frame- 
work, with fix = 1 and neglecting the radiation compo- 
nent : 



dL{z) 



du 



c(l + z) 
Ho 
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The connection with Eq. is given by : 
9W = \ + lw{z)flx{z) ; flxiz) = flx^^^^ (3) 

In Ref. 13] Riess et al. have recently presented an anal- 
ysis of 156 SNe including a few at z > 1.3 from the 
Hubble Space Telescope (HST) GOODS ACS Treasury 
survey. They use a kinematical description (no dynami- 
cal hypothesis) with a simple parametrization of q{z) in 
Eq.Qj q{z) = Qo + Qi z, and conclude to the evidence for 
present acceleration go < at 99% C.L. (go ~ —0.7) and 
for past deceleration q{z) > beyond zt = 0.46 ± 0.13. 
Concerning the dynamics in a flat universe, they con- 
clude on the validity of the Cosmic Concordance version 
of the ACDM Model that is flM « 0.3, w{z = 0) w -1 
and no rapid evolution of the EoS. 

Performing the same kind of analysis, we confirm the 
numbers and errors obtained with the Gold Sample (see 
Ref. 0) of SNe for the "kinematical fit" : qo = -0.74 ± 
0.18 and qi = 1.59 ± 0.63. In spite of the goodness of 
the fit and of the relatively small errors it yields, this 
simple strategy leads to some bias which we quantify in 
the following. 

As noticed several times, 0, IE IS H the physical 
parameters 57's and w{z) are related to the measured 
quantity dL{z) through a multiple integral relation which 



results in a somewhat uncertain determination of these 
parameters, even at their present values flM and w{Qi). 
In Ref. we have analysed quantitatively the bias which 
occurs when one tries to extract w(0) neglecting a pos- 
sible redshift dependence of w{z). In the same time, 
other fitted quantities (essentially flu) are affected, due 
to strong correlations between flM, u'(O) and dw{z)/dz. 
Another pitfall originates from the assumed value (and 
uncertainty) of the prior on flM which is used for the 
fit Q : an artificial convergence towards the Concor- 
dance Model seems to occur whereas the simulated fidu- 
cial model is very different. 

Concerning the determination of the deceleration pa- 
rameter, one encounters a similar problem : the extrac- 
tion of qo is not independent of the assumed form of the 
function q{z), a form which is related (dynamical ap- 
proach) or not (kinematical approach of Riess et al.) to 
the evolution of the DE EoS and to the value of Hm- In- 
deed, when choosing the kinematical approach, Riess et 
al. are in some sense model independent but they must 
use a particular description of the kinematics : e.g. the 
linear form q{z) = qo + qi z. It is valuable to question 
this description, for instance in the light of some dynam- 
ical models which include an evolution of the DE EoS. 

To modelize the evolution of the EoS, in the absence 
of deep physical insight, the choice of a parametriza- 
tion is arbitrary. Simplicity imposes a two-parameter 
parametrization for w{z), then consistency with the data 
from the Cosmic Microwave background (CMB) imposes 
w{z) < at high redshift. Therefore, the simple linear 
parametrization : 



w{z) — Wo + Wl z 



(4) 



we have used previously [a, |3| (with many other authors) 
is too badly behaved at high z. We prefer to switch to 
the parametrization advocated in llCl [DJ and which is 
now widely used in the literature 



w{z) = Wo + WaZ/{l + z) 

Then from Eq.©, one gets : 



(5) 



flx{z) = fix 



HHz. 



(1 + 2)^(i+"'°+"'°)e"^"'"^/(^+^) (6) 



and q{z) reads : 



1 + ^(1 + zy3(wa+w^)^3w^z/(l+z) 

fix 



(7) 



The behavior of q{z) is plotted in Fig.l for various 
models listed in Table I : it is obvious that q{z) is in gen- 
eral far from a linear shape. Even in the ACDM model, 
q{z) is not exaclty linear in z as can be seen from Eq.Q 
and Fig.l. 
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FIG. 1: q{z) from Eq.Q with Qt = 1 Eind some particular values 
of the parameters {^Im^ ^Oi ^a) given in Table I. 

TABLE I: Model examples with fiducial values {Qmt^O, Wa)- The 
fiducial go(^Mi^o) is calculated from Eq. Q. The fitted qo are ob- 
tained with a linear 3-fit {Ms, qo,qi) of the simulated data. These 
values are well in agreement with the one extracted from the Gold 
Sample (gfj = -0.74 ± 0.18). For consistency, Models A,B,C,D 
have been chosen to give a good dynamical 4-fit (Ms, '^'Oi ^a) 
of the Gold Sample. 



model 




wo 




qo{QM,WQ) 


fitted qo 


ACDM 


0.27 


-1 





-0.6 


-0.57 ±0.17 


A 


0.27 


-1 


-2 


-0.6 


-0.79 ±0.17 


B 


0.27 


-0.8 


-3.2 


-0.38 


-0.74 ±0.17 


C 


0.50 


-2.4 


1.4 


-1.3 


-0.75 ±0.18 


D 


0.50 


-0.6 


-15 


0.05 


-0.75 ±0.18 



To illustrate the consequence of the non-linear form of 
q{z) on the kinematical fit, we simulate SNe data sam- 
ples corresponding to the same statistical power as the 
true data sample by fixing the fiducial(^) values for the 
parameters Mg ,il^j,WQ and . In the following, we 
consider a flat cosmology and we fix Mg — —3.6. Then 
we perform a three-parameter (Ms, qo,qi) kinematical fit 
of the simulated data. The results for the models we use 
are presented in Table I. 

First, the ACDM Concordance model {wq — —I, 
= 0) with the value (0.27) of 0, yields a fit- 
ted value: go — —0.57 ± 0.17 which reproduces well the 
fiducial value in the limit of the errors. Note however 
that the value ^g"*" = -0.74 ± 0.18 obtained from the 
real data is la away from the actual ACDM value. 

With Model A, where a time variation of the EoS is 
allowed (w^ 7^ 0), we illustrate the fact that even if 
from Eq. Q is independent of , its fitted value is not. 

With model B (this model could correspond to one of 
the k-essence models 0|), the fitted go value is now 2a 
away from the fiducial {qg = —0.38) . 

In model C, we have changed the fiducial flf^ value to 
0.5. We have Wq < —1 (as in phantom models see 
e.g. ^] for a list of references) and a positive . The 
bias is then very large since the fitted qo is more than 3cr 
away from the fiducial value. 

One can even get a truly slowly decelerating model 
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FIG. 2: Biaised Zone (BZ) and Validity Zone (VZ) for the de- 
celeration parameter qo, for the kinematical fit {Ms,qo,qi) of 
the simulated data in the fiducial plane {wq ,w^) with Cl^ = 
0.27, Mf = -3.6. Small dark zone : where the linear approxi- 
mation q{z) = go + qiz is acceptable ; hatched zone : where the 
conclusion on the sign of go could be misleading. 

with model D : qQ — 0.05, still with fi^^ — 0.5, Wq — 
—0.6 and a very rapid and recent evolution of w{z) : 
= — 15. Of course the behavior of Model D seems 
somewhat artificial and the result of the fit can be seen 
from the non-monotonicity of q{z) in this model (see 
Fig.l), which leads to compensations in the integral of 
Eq.(^. However, let us stress again that Model D as 
well as Models A,B,C fit perfectly the present SNe data, 
yielding in particular q^ « gp"*". 

More quantitatively, it is interesting to pin down some 
regions in the parameter space where the kinematical fit 
gets in trouble. Following Ref.Ql, we fix fi^ and Mg and 
we consider the plane of the fiducial {wq,w^). In this 
plane, we define the Biased Zone (BZ) for go as the zone 
where the kinematical fit converges perfectly but where 
the difference between the fitted value and the fiducial 
value qQ is larger than the statistical error a{qo). The va- 
lidity zone (VZ) is the complementary region. It must be 
stressed that the BZ is undetectable with real data. We 
give the two zones in Fig. 2 for Of^ = 0.27. The ACDM 
falls in the VZ as expected, and the models A,B,C,D be- 
long to the BZ. Within the VZ, the zone where the linear 
approximation q{z) ^ qo + qiz of Eq.([7|l is acceptable is 
quite small. In fact the rest of the VZ is due to acciden- 
tal cancellations when Eq.© is injected in Eq.QJ. The 
lower-right corner of the BZ corresponds to models such 
that Wa ^ —SiVq^Im ; they display a non-monotonic be- 
havior for q(z) which explains the bias. In addition, even 
with = 0.27, there exists a region in this zone where 
the conclusion on the sign of qo could be misleading, i.e. 
a zone where qo + a{qo) < whereas g^ > 0. This zone 
corresponds to —0.5 < Wq < with large and negative 

Therefore, it appears that the linear expression of q{z) 
at first order : q{z) = go + gi^, cannot be satisfactorily 
used over a large range of EoS parameters. Going to a 
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TABLE II: Results of a "dynamical" 4-fit with w(z) = too + 
Waz/{1 + z) using the Gold data from Q. 



flM prior 




Wo 




*) 


0.27 ± 0.04 
0.27 ±0.20 


-0.04 

u.oo_o 21 


_1 4g+it3i 

'"-0.60 


'^•^'-'-1.60 
o rr| + 2.20 


1 19+U.33 
^•^^-0.36 

1 91+0.49 
r-^r_Q 45 



second order parametrization will not improve the issue, 
as it will only enlarge the number of parameters, will 
degrade the errors, and will not solve the bias problem. 

Alternatively, one can deduce a qo value from real data 
using Eq.(|7)) (with z = 0) where we consider the dy- 
namical description Eqs.Qand (Q, by performing a 4- 
parameter fit {Ms, nM,wo,Wa) of the Gold Sample of 0- 
We have calculated asymetric standard deviations (devi- 
ation Sx^ = 1 for a given parameter), marginalizing over 
the other parameters by minimisation. For go: which is 
a derived parameter, we use a Monte-Carlo technique : 
generating 1000 simulated experiments around the data 
best fit, we estimate the width of the go distribution. 

The extracted values are displayed in Table II. In this 
Table, we have used a strong prior on Hm (0.27 ± 0.04) 
and also a more reasonable weaker uncertainty of ±0.2, 
since adopting strong priors is recognized to be dangerous 
(see 1^ U^)- Comparing with our previous fits, we see 
that changing the w{z) parametrization from Eq.Q to 
Eq.ljHJ has not changed the behaviors discussed in 
In particular choosing a strong prior around = 0.3 
pushes wq towards -1 whereas, with a weaker prior the 
central value of wq is clearly lower than -1, a fact which 
has been noticed by many authors (see e.g. 0, 0, fSH 

El 113). 

The estimated go values confirm that the Universe is 
presently in an acceleration phase, at 95% C.L. We note 
that the extracted value for a weak prior is 2a away from 
the ACDM value go = —0.6. This is a direct consequence 
of the extracted value of the point {wo,Wa), which is 2a 
away from (—1,0). With the strong prior, the errors on 
go, as on other parameters, are systematically smaller 
when is close to the Concordant Model value (w 
0.3). This conclusion would be different if we use a larger 
central value for ^Im ■ in this case, the go value would be 
even more negative, but with such a large error that it is 
not possible to conclude. 

Finally, if no prior is applied on Qm , one can conclude 
with some caution that go < at 80% C.L. 

We can also address the problem of the determina- 
tion of the value of the transition redshift zt between 
deceleration and acceleration. Since the kinematical fit 
yielding the go and gi values found in Ref.0| is biased, 
one cannot be too confident with the advocated value 
Zt = 0.46 ±0.13. 

Fig. 3 displays the function g(z) — qbf{z) computed 
from the best fit {bf) values of the parameters enter- 
ing Eq.0. The la and 2a intervals are computed by 
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FIG. 3: best fit qi,f{z) function computed with Eq.JTJ using the 
best fit values of the cosmological parameters from the current data 
(heavy plain curve). The uncertainty contours are constructed from 
the probability density function around the best fit at a given red- 
shift from one thousand simulated experiments. The dark stripe 
corresponds to the Icr error on the fitted function from a simula- 
tion around the ACDM concordant model fiducial values, with the 
SNAP statistics. 

taking the probability density of simulated experiments 
around the qbf{z) value in each redshift bin. We still 
use a weak prior Hm = 0.27 ± 0.2. The extracted Zt 
value corresponds to the point where qhf(zt) = and 

%^(z = Zt) > 0. We find zt = 0.34+^' J^, where the er- 
rors are evaluated from the Zt probability density. With 
the dynamical g(z), our errors are small because the 
transition is steeper than in the linear kinematical case 
(g(z) = go ± giz), as can be seen from Fig. 3. 

Our central value is far from the " theoretical" Zt of the 

Concordance Model : Zt = [^V^^ - 1 = 0.76 (for^M = 

0.27). Moreover, we have observed that, contrary to the 
extraction of go, the extraction of zt is very much depen- 
dent upon the chosen parametrization for w{z), a fact 
mentioned by various authors [23, ll^l- Then, one can 
stay prudent and state that, with present data, the ex- 
tracted Zt value is certainly much less robust than the 
extracted go value, even in a dynamical approach. We 
should also point out that no systematical effects are 
taken into account in these evaluations, in particular a 
normalisation effect between low and high redshift super- 
novae would affect directly the determination of Zf. 

It is interesting to evaluate what should be the best 
strategy to extract go and Zt in the future. Using a bi- 
ased method, as the linear kinematical one, will be worse 
when the statistical errors will decrease. A dynamical, al- 
though model-dependent approach, remains a good way 
to estimate the behavior of the acceleration. In Fig. 3, 
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we show the results for a statistical sample expected from 
the space-based SNAP mission where the systemat- 
ical uncertainties should be controlled at the same level 
of precision. Thanks to the rich sample of 2000 SNe and 
also to the large range of redshift (0.2 ^ z < 1.7) a very 
precise determination of the transition region should be 
allowed. For instance, with the ACDM fiducial values 
(wq — —l,Wa — 0,Q,M = 0.27), keeping a weak prior on 
Hm, one gets Zt = 0.67t^:EJ^ and qo = -O.bbtolt- 

We estimate that the correct procedure for the future 
is to avoid using priors on fijvf- We prefer instead the 



combination of supernovae data with other probes such 
as the CMB, the Weak Gravitational Lensing and the 
galaxy power spectrum. This kind of method will not 
introduce external biases and sonic recent papers go in 
that direction (see e.g. ,16, ■21.. .22. .25^). 
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